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Abstract 

Motivated by recent discussions concerning the manner in which damping appears in the 
electric polar izability, we show that (a) there is a dependence of the nonresonant contri- 
bution on the damping and that (b) the damping enters according to the "opposite sign 
prescription." We also discuss the related question of how the damping rates in the polar- 
izability are related to energy-level decay rates. 
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I. INTRODUCTION 



Several recent papers address the question of how material damping effects should 
be included in the response of an atom or molecule to an applied electric field 
We will consider the simplest case, that of the linear atomic polarizability, which in 
the absence of damping is given by the Kramers-Heisenberg formula, 

a,{uj) = 4 y knl' i + (1) 

for state i. Here uji and rji are the transition (angular) frequency and coordinate 
matrix element, respectively, between states i and j, and the field frequency u is 
assumed to be far removed from any of the atomic transition frequencies uji. More 
generally one associates damping rates 7jj with the different transitions and writes 

a,{uj) = 4 y knl' ( + ^—^^) (2) 

where ^ = +1 according to the so-called "opposite sign" prescription and ^ = — 1 
in the "constant sign" prescription. The difference appears only in the nonresonant 
terms, and is therefore unimportant in most situations. However, the question of 
which prescription is the correct one raises some interesting points, as we shall see, 
and the purpose of this paper is to address some of these points as well as to answer 
the question of whether one should take ^ = -|-lor^ = — lin equation (jj)). 

One might ask first whether a damping term should appear at all in the nonres- 
onant part of the Kramers-Heisenberg formula, i.e., whether we should in fact take 
^ = instead of either ^ = +1 or ^ = — 1. An analysis involving the diagonalization 
of the (two-level) atom-field Hamiltonian in the rotating-wave approximation, for in- 
stance, shows that there is no damping term in the nonresonant denominator ^, a 
result that is certainly accurate for most practical purposes. In a broader context the 
issue here is an old one. Thus the imaginary part of the polarizability Q implies an 
absorption coefficient having the usual Lorentzian form 



as well nonresonant part 

(4) 

for a transition of frequency and linewidth 7, and one might question whether, as a 
matter of principle, Q contributes to the absorption lineshape. In his consideration 

n 

of possible corrections to the Weisskopf-Wigner lineshape. Lamb [7] noted that "such 
a contribution [as (jH)] appears in some derivations," but added that it would be 
negligible compared with the resonant contribution Q. 

The effect of damping on the nonresonant part of the polarizability is not an 
entirely trivial matter, and the literature relating to the subject reveals significant 
disagreement on some rather basic aspects of dissipation theory. The purpose of this 
paper is to address the principal points where there is disagreement and to obtain 
what we regard as the correct form of the polarizability when damping is included. 

In the following section we consider the problem of the electric-dipole interaction of 
a two-level atom with the quantized electromagnetic field, assuming that all but one 
of the field modes are initially unoccupied. Using the rotating-wave approximation 
(RWA) for the atomic source field but not for the applied field, we obtain exactly 
.he ..ult c.ed ea.ie. 0, and i„ particular we End .Ka. the.e is no contribution 
from (radiative) damping to the nonresonant term in the polarizability. In Section 
mil we go beyond the RWA in the atomic source field and find that the damping now 
appears in the nonresonant term, and that it does so in accordance with the "opposite 
sign" prescription. Section HXl presents a discussion of these results, including their 
connection to the classical theory of radiative damping. Section |3 focuses on the 
form of the damping rate 7jj, and we argue that, contrary to what sometimes appears 
in the literature, 7jj depends on the sum rather than the difference of the level decay 
rates. Our conclusions are summarized in Section IVll 
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II. DERIVATION OF LINEAR POLARIZABILITY: RWA 



The model we consider is described by the Hamiltonian 

H = hjjQ(Tz + ^ fkUkolak - dEa^ = hcooCTz + ^ hjJka\ak - ih ^ - «!) (o" + cr^) 

k k k 

(5) 

for a two-level atom (TLA) with transition frequency ujq and dipole moment d inter- 
acting with the electromagnetic field. The a operators are the usual Pauli operators, 
with a and o"^ the lowering and raising operators for the TLA. and a|. are the anni- 
hilation and creation operators for field mode k, and = (d ■ A//i)(27r/lco'fc/l^)^/^, 
with V the quantization volume. The subscript k denotes (k, A), where k is the 
wave vector assocated with a plane- wave mode of frequency u^. = |k|c and e^^A is a 
corresponding polarization unit vector (k ■ ek,A = 0, ek,A • a' — '^aa', A = 1, 2). 

The commutation relations for the atom and field operators give the Heisenberg 
equations of motion 

a = -iuoa + ^ Ck{cr^ak - a^cr^) (6) 



-iujkak + Ckicr + (t"^) (7) 



We have chosen a normal ordering for the field annihilation and creation operators, 
which is especially useful in the case that the applied field is described by a coherent 
state [Eq. (fTT| ]. As we are interested only in determining the linear response, the 
equation of motion for will not be needed for our purposes. 
The formal solution of equation ((7j) is 

(t) = afc(0)e-'"'=* + Ck f dt'[a{t') + at(t')]e^"^ (8) 



afcl 

In one version of the rotating-wave approximation (RWA) we ignore the coupling 
between the creation operator for the field and the raising operator for the atom; 
this corresponds, for the purpose of obtaining the equation of motion for the field 
operators, to the neglect of the terms a\a'^ and aa^ in the Hamiltonian (jSj). In this 
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approximation the equation of motion for a becomes 

k 

+ Cl f dt'a,{t)a{ty^>''^''-'^ (9) 

Note that we are not making an RWA in the free-field operators 0^(0) and a|.(0), so 
that both annihilation and creation free- field operators [afc(O) and a|.(0)] appear in 
©• 

Taking expectation values over the initial atom- field state on both sides of 0, we 
have 



k 



k -^0 



(10) 



We assume that the initial field state \iPf) corresponds to a single occupied mode 
described by a coherent state with 

afc(O)lV'F) = a\ijjF) , (^^^14(0) = Oi*{i'F\ (11) 
corresponding to the expectation value 

/ 9 ft \ 1/2 

{E{t)) = I [ j [ae"*"* - a'e'^'] = So cos cut (12) 



of the applied electric field. Thus 



5^C,[(a,(t)afc(0))e-*-'^* - (4(0)a,(t))e*-'^i = - -d^S^ cos ut {a, {t)) (13) 
k 

where dx is the component of the dipole matrix element along the direction of the 
applied field; rf^ = d'^/3 for the spherically symmetric atom. 

We assume that the operator a^, which corresponds to the population inversion, 
changes sufficiently slowly that we may take 

(a.(t)a(tO) = {cr.{t')a{t')) = -{a{t')) (14) 
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in the integral appearing in (jlUj) . Since we want to obtain the polarizabihty for 
the TLA in a particular state, we assume further that the atom remains with high 
probability in its initial state. Assuming this initial state to be the lower state, we 
approximate {az{t)) by —1, so that, using ()13|) and the approximation ()14j). we replace 

m by 

^ ft 

{&{t)) = -tuo{a{t)) + i-d,Socosut -Y.Cl / dt\a{t'))e"^^^''-''> (15) 

We seek a solution of ()15|) of the form 

{a{t)) = se-*"* + re''^' (16) 

with s and r constants to be determined. This implies 

i{uJo - uj)se-"^' + i{uJo + uj)re'^' = -^So{e-'^' + e*"*) - [7_(c^) - iA_(cu)]se-*"* 

-[7+(o.)-^A+(o.)]re-* (17) 

where 



u Jo 67r nv J ^ 

de sin^ e 



00 



k 

d' 

'2d'' r , . 2d2o;3 



and 



Jo 3n.c 



for t >> l/ciJ, where U is the unit step function. Note that the damping rate 7-(co') is 
frequency-dependent (A-|-(ti;) is obviously divergent but, as discussed in Section 
4, this has no direct bearing on our conclusions regarding the effect of damping on 
the polarizabihty.) To obtain the polarizabihty a{uj) we write 

P = d,{a,) = d^{{a) + (at)) = 24■7^e[(r + s*)e-^"*] = ne[a{uj)Soe"'''] (20) 
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for the induced dipole moment. This yields 

^ f 1 ^ 1 

3^ — lu — A_(ct;) — i7-(co') uq + uj — A+(co') + «7+(co') 

Note that 7+(co') = 0, and that therefore there is no damping contribution to the 
second (nonresonant) term. 7-(to'o) is half the radiative decay rate of the upper state 
in the absence of any applied field. 




III. DERIVATION OF LINEAR POLARIZABILITY WITHOUT RWA 

Let us now recalculate the polarizability, this time retaining both terms inside the 
integral of equation (jSJ, i.e., without making the RWA in the (source) field produced 
by the atom under consideration. Then (jH)) is replaced by 

(6"(t)) = — iu^ia it)) — i^8o cos ujtiozit)) 

n 

+ E^' f dt'[{az{t)a{t')) + {a.{t)a\t'))Y^^^'-'^ 

k "^0 

- E f dt'[{a\t')a-M) + (a(t')^.(t))]e-^-(*'-*) (22) 

when we take expectation values as before. The approximations tantamount to (fT^ 
are 

(a.(t)a(O) = {az{f)a{f)) = -{a{t')) 
{a^a^t')) - (a.(t')at(O) = {<r\t')) 

{a{t')a-M = Wit>.{t')) = Wit')) (23) 

where we use the equal-time identities crz{t)(T(t) = —a(t)(Tz(t) = —<j(t). Using these 
approximations in (j22I), together with the approximation (cXzit)) = — 1 in the second 
term, we obtain the non-RWA extension of 

i&it)) = -iuj^{a{t)) + '-fS^cosut + y^Cl / dt'[-{a{t')) + {a\t'))]e'^^^''-'~^ 

-Y,Cl f dt'[-{a\t')) + (a(t'))]e-^-'=(*'-*) (24) 

^0 



It is important to note that in equations (j23|) we have used the commutation relations 
between az{t) and o"(t), and have obviously not made the approximation that 

cr^ could be replaced by —1. The latter approximation is made only in the second 
term of ()22|) . where cr^ multiplies the applied field but no atom operator, so that the 
approximation does not violate the commutation relations from which we obtained 
the equations of motion. The two approximations are different: that made in (j23|) 
assumes that ^^(t) varies little on time scales ~ l/cj^ for field frequencies uJk ^ uj that 
will contribute significantly to the variation of {a{t)), whereas that made in replacing 
{az{t)SQ cosut by —Sq cosut assumes that the atom remains with high probability in 
its lower state because the field frequency lies outside the absorption linewidth. The 
difference between these two approximations involving az turns out to be irrelevant 
for the final results when the RWA is made, as is clear from ()14|1 . 

We again have a solution of the form now with s and r satisfying 

Xs + Ur* = 

Vs + Yr* = ^So (25) 

where 

X = ujo-u- [A_(cu) - A+{uj)] - i[7-(cu) + 7+(cj)] 
U = [A_M-A+H]+*[7_M+7+M] 

Y = UJO + UJ+ [A^iuj) - A+{uj)] + i[-f^iuj) + 7+(^)] 

V = [A4uj) - A^Lu)] + ^[7-(^) + 7+(^)] (26) 

Assuming that 'J±{uj) and A±{uj) are small in magnitude compared to cjq i u, we 
have 

dxSo 1 dxSo 1 



2h X 2h ujQ-uj- \A_{uj) - A+(u;)] - i[i-{uj) + l+{uj)] 

dx^o 1 _ dxSp 1 

~Wy ~ 2h iUo + iu+ [A_(cu) - A+(cj)] + t[^-iuj) + 7+(^)] 
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and, from (pUj) . 

a{uj) - ^^^^^ 

uo-u - [A_(u;) - A+{u)] - i[7-(u;) + 7+(^^)] 



UJO + UJ+ [A_(cu) - A+{uj)] + i[7-(cu) + 7+(cj)] 



(2^ 



IV. DISCUSSION 



In contrast to the RWA result (|2ip. A±{uj) and 7±(c<j) appear in both the resonant 
and nonresonant terms of (j28j) . Consider first the physical significance of A-|-(ci;), 
assuming that the frequency u of the initially occupied field mode is sufficiently close 
to ujq that we may take 

and focusing only on the resonant term in aiiS). In a more complete analysis involving 
the transformation from the fundamental minimal coupling form of the Hamiltonian 
to the electric dipole form, it is found that the additional term 27r f d^P-*- • P-*- ap- 
pearing in the transformed Hamiltonian has the effect of replacing ()29|1 by j2| 

A4^o ^ / 30 

With this modification it is seen that A(co'o) = A_(co'o) — A+(co'o) is simply the (un- 
renormalized) TLA radiative frequency shift, i.e., the difference in the radiative level 
shifts of the two levels 1^ . In general, however, the approximation (j29p is not appli- 
cable, and the radiative level shifts 1iI\±[lS) depend on the frequency of the initially 
occupied mode. In the polarizability (PHj) the frequency shift A(co') = A„(co') — A+(co') 
adds to the field frequency lo in both the resonant and nonresonant terms, whereas 
in the RWA A^iuS) does not appear in the resonant term and A_(co') does not appear 
in the nonresonant term. In other words, the RWA does not correctly include the 
radiative frequency shift as the difference in the radiative level shifts of the TLA. 

The expressions for the level shifts hI\±[iS) are specific to the TLA model, but 
are easily generalized to the case of a real atom. This extension, even with the 
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standard renormalization procedures, still leaves us with divergent level shifts in the 
nonrelativistic approximation. A high-frequency cutoff mc^/h results in Bethe's ap- 
proximation to the Lamb shift Since this procedure is very well known, and 
we are in any case only concerned with the form in which the radiative corrections 
appear in the polarizability, and not their numerical values, we will simply assume 
henceforth that the frequency shift has been accounted for in writing uq ± uj. 
Thus 

, , cP f 1 1 \ 

a(uj) = — + 31 

3n \ujo — UJ — 27(0;) ujq + UJ + i'~i\uj) J 

where "^{uj) = 'y-{uj) + 7+(co'). Like A(u;), ^7(^0') is effectively an addition to the 
applied field frequency uj. Unlike the frequency shift, however, the damping rate '~f{uj) 
is half the sum of the decay rates 'y±{uj) of the two levels. Of course the decay rate 
7+(c<j) of the ground state in our two-level model is zero but, as discussed in the next 
section, (j3ip is valid more generally when the decay rate of the lower level of the 
transition is not zero. That is, the damping rate appearing in the contribution to the 
polarizability from any given transition involves half the sum of the decay rates of 
the two levels of the transition. 

Regardless of whether the lower-level decay rate vanishes, the non-RWA result 1)311) 
shows that both the resonant and nonresonant contributions to the polarizability have 
a nonvanishing damping term in their denominators, this damping term being half the 
upper-level decay rate. In particular, it is seen that the damping appears according 
to the "opposite sign prescription," i.e., ^ = +1 is the correct choice in the dispersion 
formula The same conclusion was reached by different lines of reasoning by 
Buckingham and Fischer Q]. 

Note that, if 7 is taken to be a (positive) constant, independent of frequency, then 
the opposite sign prescription is consistent with the causality requirement that the 
polarizability should be analytic in the upper half of the complex uj plane [ll|. But in 
general the decay rates are in fact frequency-dependent ^, and causality is ensured 
only if the model used to calculate ■y{uj) is itself causal. In fact, as recalled below, 
radiative damping provides an example in which this is not the case. 
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In one approach to a classical calculation of the natural lineshape, one considers the 
solution x(t) = Aqc"'*'* sin(ci;ot) of a damped dipole oscillator with resonance frequency 
ujq. The lineshape is taken to be proportional to the squared modulus of the Fourier 
transform 

a{ij) = ^ [ rfte-^*e'"*sin(cJoi) oc ( + ) (32) 



27r Jq \ujq — uj — i'f ujq + uj + i'y 

and is seen to be consistent with the "opposite sign" prescription. In contrast to this, 
an old paper by Weisskopf ^ | implies the result 

a(a;) oc ( ^ (33) 

\tUo — LU — i'j coq + u — i'j J 

which is consistent with the "constant sign" prescription. However, since this result 
is based on the integral appearing in ()32|). it seems that involves a sign error or 
perhaps just a typographical error. 

Since the absorption coefficient may for our purposes be taken to be proportional 
to the imaginary part of a{uj), equation (p?T|) implies an absorption lineshape propor- 
tional to 

^(^) = 7 ^^\2^ — 2 ~ 7 1 \ I — 2 (34) 

{ujo - Luy + Y {^0 + t^) + 7 

The same result, for 7 taken to be a constant, was obtained on the basis of the 
Lorentz model by Van Vleck and Weisskopf who noted that the minus sign in 
the nonresonant term "must be used because the excitation of the molecule is here 
accompanied by emission rather than absorption of a light quantum," a process which 
is excluded when the RWA is made |l4 1 . 

It is also of interest to compare the result (jHTj) with the corresponding result given 
by the classical theory of radiative damping based on the equation 

X + UqX — -x=e£ocosujt (35) 



3mc^ 

The polarizability of the classical dipole oscillator described by this equation is 

e^/m 

Oid[^) = —5 ; 2/- 9 / vT^ 

uJq — uj'^ — ^[ze^/mc^juj-^ 

/ 1 1 A 

+ , . , ... . .^ (36) 



2mujQ \uJq — LJ — ijcii^) uJq + LU + i'jcii^^ 
11 



where lu'q = a/cUq ^ 7cz(^) Idi^) = (e^/Smc^)^^^. The replacements /2mujQ 
e^/i/2ma;Q and e^oj"^ /?>mc' — > 6^/2^0'^ /3mc'^, where /i = 2mujQ(f/e^h and /2(i^) = 
{2muj(P / e^h), make the classical result ()35p equivalent to ()31|). These replacements 
involving effective oscillator strengths /i and /2 are the usual substitutions required to 
put classical oscillator results in agreement with some of the corresponding quantum- 
mechanical expressions. 

The uj'^ in the denominator of (|36|). or in other words the third derivative of x in 
equation (j^ . leads to a pole in the upper half of the complex oj plane, thus violating 
the causality requirement that the polarizability be analytic in the upper half-plane. 
The nonrelativistic theory of radiative reaction is well known to be acausal, but the 
acausality occurs on such a short time scale that relativistic quantum effects must 
be taken into account. For most practical purposes the acausality is of no conse- 
quence. Thus, for instance, equation (j36j) leads to the correct extinction coefficient 
(oc ti;Tm[Q;(ti;)]) due to Rayleigh scattering. 



V. RELATION OF DAMPING IN THE POLARIZABILITY TO LEVEL DE- 
CAY RATES 

These considerations are easily extended beyond the two-level model, with the 
result that the linear atomic polarizability has the form 

«^(^) = 4 E i^^-^i' { — + ^ i , (37) 

The damping rate 7ji(co') has a "dephasing" contribution associated, for instance, 
with elastic collisions, as well as a contribution associated with the decay rate of the 
atomic states i and j. Here we consider only the latter contribution, which is due to 
radiative decay and other loss processes. In the case of radiative decay, for instance, 
7jj is found, by a straightforward multilevel generalization of the calculations in the 
preceding sections, to be half the sum of the radiative decay rates associated with the 
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two states i and j: 




Ei>Err 

i3 Vitr 



where denotes the energy of state m. If we replace u by lJq in this formula, we 
obtain half the spontaneous decay rate of state j in the case that the field is initially 
in the vacuum state. This result was obtained, for example, by Weisskopf and Wigner 
jlS^ . Landau Q|, and many others 0]. The same conclusion is reached in the more 
general case where the energy levels decay by nonradiative channels: 'yji{uj) is half 
the sum of the total decay rates of the states i and j. 

Various authors, however, have calculated or assumed — erroneously, in our 
opinion — that 7jj involves the difference in the decay rates of the states i and j [Ij 
jisl ] , [l^ . In addition to the Heisenberg-picture calculation leading to the conclusion 
that 7jj involves the sum rather than the difference of energy-level decay rates, as 
presented in this paper, the following simple argument can be used. Let Cj(t) and Cj{t) 
be the (Schrodinger-picture) probability amplitudes for states i and j, and let 7^ and 
7j be the decay rates of these states. Then c*{t)cj(t) and Ci{t)c*{t), which determine 
the polarizability, decay at exp[— |(7j+7j)t], and so the linewidth in the polarizability 
must involve the sum of ■ji and 7j rather than the difference. Sushchinskii ^|, for 
instance, expresses his results in terms of complex energies E'- = Ei — \iTi and their 
differences E[ — Ej = Ei — Ej — |i(r,j — Tj), whereas the appropriate differences 
entering into the polarizability are E^* — E'- and E[ — E'*. 

Finally we note that Andrews et al. [19] have stated a polarizability sum rule which 
in the simplest case of the linear polarizability can be expressed as ai{uj) = 0. A 
physical plausibility argument for this sum rule can be adduced as follows. If pi is the 
probability that the atom is in state i, then the linear polarizability at field frequency 
uj is 

a(t^) = ^Piai{uj) (39) 

i 

Consider the idealized limit in which all the Pi are equal. Then the polarizability and 
therefore the induced emission or absorption rate at frequency uj becomes proportional 
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to just J2i Q^j(^)- But if all the states are equally populated the net induced emission 
and absorption rate must vanish, implying the polarizability sum rule conjectured by 
Andrews et al. From the expression (jHTj) it follows that this sum rule is statisfied 
only if is symmetric in i and j, i.e., must involve the sum rather than the 
difference of 7^ and 7^. (We note that, in the case of the constant sign prescription for 
the damping terms in the polarizability, the polarizability sum rule would be satisfied 
only if 7ji were antisymmetric in i and j.) 



VI. SUMMARY 



Following a standard, nonrelativisitic approach, we have considered specifically 
the case of a two-level atom interacting with the quantized electromagnetic field, one 
mode of which is initially occupied and described by a coherent state. Working in 
the Heisenberg picture, we calculated the polarizability with and without making 
the RWA for the atomic source field. In the RWA we obtained a known result, 
and in particular the nonresonant contribution to the polarizability was found to 
have no damping factor in its denominator. Going beyond the RWA, however, we 
found that both the resonant and nonresonant contributions to the polarizability have 
the radiative damping rate in their denominators, and that the polarizability has a 
form that is consistent with the so-called opposite sign prescription for including the 
damping. 

The radiative frequency shift appearing in the non-RWA expression for the 
polarizability depends on the radiative level shifts in the correct way, i.e., it is the 
difference of the two level shifts. The damping rate appearing in the non-RWA 
expression for the polarizability is half the sum of the radiative decay rates of the 
two levels, in contrast to the difference of the decay rates that has been obtained or 
assumed in some treatments. The fact that the polarizability depends symmetrically 
on the decay rates of the energy levels is consistent with the polarizability sum rule 
of Andrews et al. [3] when the (correct) opposite sign prescription is used. 
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